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Abstract 

We compute as a function of rapidity y the baryon number carried by quarks and 
antiquarks with p T > p ps 2 GeV produced in Pb+Pb collisions at TeV energies. The 
computation is carried out in lowest order QCD perturbation theory using structure 
functions compatible with HERA results. At po = 2 GeV the initial gluon density is 
both transversally saturated and thermalised in the sense that the energy/gluon equals 



X 



to that of an ideal gas with the same energy density. Even at these high energies the 



initial net baryon number density at y = at r = 0.1 fm will be more than the normal 
nuclear matter density but the baryon-to-entropy ratio is only (B — B)/S ~ 1/5000. 
Further evolution of the system is discussed and the final baryon-to-entropy ratio is 
estimated. 
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1 Introduction 



The formation of QCD plasma in central collisions of large nuclei in the TeV energy 
range can be described as the formation of a large number of partons (gluons and 
quarks) with transverse momentum p^ > Po, where po parametrises the lower limit of 
the pt of the partons included in the computation [|J . The number of these partons 
in an average collision can be estimated in QCD perturbation theory H,f|. It scales 
~ A 4 / 3 and for p « 2 GeV its value is easily (9(1000) per unit of rapidity in the central 
rapidity region. This process takes place in a time l/p ~ 0.1 fm/c and it is thus equal 
to the initial formation of QCD plasma. Making predictions for the final observables 
requires modelling the further evolution of the system either hydrodynamically j|, |(| 
or in various kinetic models [|7||-||13||. 

For a thermal system the computation of the overall multiplicity of the initial partons 
gives the entropy, associated with the temperature T(t,x), and a computation of the 
quark and antiquark number gives the net quark number, associated with the chemical 
potential /z(t,x). The purpose of this article is to compute in more detail the initial 
amounts of gluons, quarks and antiquarks with p^ > p Q and from this to estimate the 
final baryon-to-entropy ratio. 

At LHC energies, \fs= 5500 GeV, the gluons are by far the dominant component 
of the initial state. In fact, for po = 2 GeV we shall see that QCD perturbation 
theory implies that the gluonic subsystem has two remarkable properties: it is both 
thermalised, in the sense that its energy/gluon is the same as in a thermally equilibrated 
ideal gas with the corresponding energy density, and transversally saturated |TJ, [I], [J] in 
the sense that the total transverse area occupied by the gluons equals the total nuclear 
area. Thus, although p in principle is a parameter constrained by soft physics JL6|, [17| , 
its value 2 GeV gets a well defined dynamical significance in heavy ion collisions at 
very high energies. 

After the initial production the system goes through various expansion and phase 
transition stages. A natural first assumption is to assume that the expansion is adia- 
batic, so that both the total entropy 5* and the net baryon number B — B in a unit 
of rapidity are constant. This, in particular, implies that Et in a unit of rapidity de- 
creases rapidly, the work done against expansion transfers Et to larger rapidities from 
the central unit. We shall analyse this decrease, as well as possible increase of S and 
B — B due to dissipation and non-perturbative source terms in some detail. 

Concerning B — B, the system produced in pp collisions is different from the one in 
AA collisions. In pp the color correlations (strings) between individual quarks (anti- 
quarks) produced in the central rapidity and the (di-)quarks left in the fragmentation 
regions are essential for baryon production [12|, [13], || |l(J. The system is dilute and 
essentially free streaming, and the (net) baryon number remains mostly in the fragmen- 
tation regions, while the mid-rapidity quark(antiquark) usually becomes an ingredient 
of a meson. Therefore, in pp collisions the initially produced q — q distribution in the 
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mid-rapidity range does not directly reflect the final net baryon distribution. In heavy 
ion collisions, however, the situation should be different. At very high energies, the pro- 
duced parton system in the central rapidity region is extremely dense, and its further 
evolution depends rather on the secondary collisions of partons (and strings as well) 
than on the direct color correlations with the fragmentation regions. This dynamics is 
built in hydrodynamical expansion equations, which imply in the adiabatic case that 
B — B in a unit of rapidity is constant. 

At \fs= 5500 GeV the total number of quarks and antiquarks will be seen to be 
initially only 10% of that of gluons. The net quark number is much less but still 
calculable. In fact, we estimate that the initial (B — B)/S ~ 1/5000 at r = 0.1 fm 
which then may increase by up to a factor 4, while in the a/5=20 GeV range numbers 



of the order of 1/20... 1/50 have been observed |]T8 , 19fl . Our result is in qualitative 
agreement with the expectation that at very high energies the baryon number of the 
initial nuclei goes to large forward and backward rapidities leaving the central region 



essentially depleted of net baryon number pi EIL Experimental data at AGS (v / s=5 



GeV) and SPS (i/i = 17 GeV) |E| have shown that nuclei are much more efficient in 
stopping baryons than pp collisions and it has been recently argued in a gluon field 
model of baryons p3| that this would persist to LHC energies (a/s in the TeV range). 
Also effects on the baryon-to-entropy ratio due to event-by-event fluctuations of the 
net baryon number have been discussed recently |^4| . 

This approach of computing gluon and quark production stretches QCD perturba- 
tion to its limits and the sources of error are obvious: next-to-leading order (NLO) 
corrections with the present acceptance criteria have not been computed, there are un- 
certainties associated with small-x gluon shadowing and with modeling nuclear effects 
for gluons in general, the quantitative results will somewhat depend on the parameter 
Po defining the lower limit of the pt of quarks and gluons included. We offer the results 
for what they are worth; they will anyway set the stage for further refinements. 



2 The formulas 

The lowest order (LO) formula for inclusive production of partons in pp collisions can 
be found e.g. in Refs. ||25|| . Our main focus is on the production of semihard partons, 
defined by pr > Po- The cross section for each flavour / can be written as follows: 
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(Vs,Po) = J dp\dy 2 Y^ x ifi/p( x uQ 2 ) x 2fj/p(x 2 ,Q 2 
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where the fractional momenta of incoming partons are 



Xl = ^(e y + e m ) and x 2 = ^(e~ y + e~ V2 ), (2) 
' s \ s 



the region of integration is 



P 2 <P 2 T <— 4l- and In^-e^^y^ln^-e^), (3) 
4 cosh y Pt Pt 



where y is fixed so that 



M < lr /v^+ Is _ \ 
m ~ \2p y 4p§ J y) 

and the invariants at the parton level are 

t = -p 2 T [l + e -fe-^)] and u = -p%[\ + e ( ^ 2) ]. (5) 
The sum over initial states includes all combinations of two parton species: 

u = gg, gq, qg, gq, qg, qq, qq, qq, qq, (6) 

while the final states consist of all pairs (without a mutual exchange): 

(kl) = gg, gq, gq, qq, qq, qq- (7) 

The factor 1/(1 + Ski) is a statistical factor for identical particles in the final state. 
In the initial state, three (anti)quark flavours are taken into accout: q = u,d, s, and, 
in the final state four (massless) flavours: q = u,d,s,c. In the running a; s (pr) we ? 
accordingly, use Nf = 4. Since 

~dt^ ,u)= £ M () 

there are only eight different types of subprocess cross sections: 

gg -> gg gq^ gq 
gg -> qq qtqj -> q%qj, i ^ j 
qiqi -> ftft ft ft -> g^g?, i ^ j 

ftft -> ftft ftft -> S-flS 



which can be found, e.g., in refs. 

For the parton densities fj/ p , we use the set GRV 94 LO |26| with a scale choice 
Q = Pt- With this set we can consistently extend our lowest order computation even 
below px = 2 GeV. 
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Summation over final state parton flavours g, u, d, s, c now gives the production cross 
sections da 9 /dy, da q /dy, da^/dy of gluons, quarks and antiquarks with px > Po in a pp 
collision. These are normalised so that 

dV lLy " 2<7 ^d(Vs,Po) (9) 
and 

°"hard(V^)Po) + 0"Ld(V^)Po) + ^hard^Po) = ^hard ( V 7 ^, ) , (10) 

where Oha^V^Po) is the cross section for one (semi)hard collision and the factor 2 
comes from two partons in one (semi)hard event. Since we will be especially interested 
in the central rapidity region, \y\ < 0.5, let us also define for later use 

1 /•+ - 5 d&f 

0hard(V*:Po, \v\ < 0-5) = - / dy — (11) 

2 J-o.5 dy 

The next-to-leading (NLO) order corrections to eq. QTJ) can be large and UAl minijet 
data indicates that a K factor of about 2.5 is needed at those energies. NLO jet 
analysis shows |27j that at higher energies the factor is somewhat less, K« 1.5. To get 
a conservative estimate we shall use K = 1. 

After the cross sections are computed on the pp level, they are transformed to the 
number N of events of the type described by a in an average central (zero impact 
parameter) A+A collision by multiplying by the nuclear overlap function TAA(b) at 
b = 0: N = crTAA(b = 0). For example, the number distributions of produced flavour 
/ are obtained as 

dN f „ .do 1 , , 

_ = TAx(b = 0) _ (12) 

and the total number of partons of flavor / within \y\ < 0.5 as 

C r hard(\/s,P0, \v\ < 0.5) (13) 

Numerically Taa(0) ~ A 2 /nR A and we shall use the value 

32 

Tp b p b (0) = — . (14) 

This procedure is equivalent to neglecting nuclear shadowing, i.e. assuming that 
fi/A(x) — Afi/ P (x). When estimating various quantities we shall further use 

.Rpb = 6.54 fm, ixRl h = 134 fm 2 , V t = 7rRl h Ay/p = 13.4 fm 3 , (15) 

where Vi is the initial volume at the time T{ = 1/po (with p = 2 GeV) when the quarks 
and gluons were formed. Note that for the longitudinal width for the volume near the 
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central rapidity, we have approximated Az »s TjAy, i.e. we have neglected effects of 

the longitudinal spatial spread of the beam-partons. 

A lower limit for the parameter p Q is obtained from the transverse saturation criterion 
|TJ|, |J. The partonic subsystem is saturated if the total transverse area occupied by 

the produced partons = NaaO° = 0)n/p 2 ] is larger than the total nuclear area = irR A - 

For Pb+Pb and iVpbPb(b = 0) = 2cJhard(v / ^;Po; \y\ < 0.5)T Pb pb(0) this converts to 

rnb > 67 v2GevJ " (16) 

The formulas required for the computation of .Ex-distribution in a hard collision 
with an acceptance region can be found in 0. We define our acceptance through a 
"measurement function" e(y), which will in the following be chosen as a step function 

, v = f 1, if \y\ < 0.5 (u) 
e \y) _ I 0, otherwise. ^ ' 

In LO the (mini)jets are produced back-to-back in px, and the -Bp-distribution in the 
chosen acceptance window can now be defined for each flavor / as 



da? If 1 
-TEr = „ / dp%dy x dy 2 si/i/ P (si, Q 2 ) x 2 f j/p (x 2 , Q 2 ) x 



(ki) 



with the normalisation 



(da ij ^ kl 

x l ~dl ^^"> 6 ^ Et ~ i 6 fk e (yi) + VMW) + 

+-g (u,t) 5(E T - [S fl e( yi ) + 5 fk e(y 2 )]p T )j, (18) 

/do 
dET dE~ = ahaTd (^> p °)- ( 19 ) 



As explained in 0, this formulation takes into account that of the two minijets pro- 
duced in the hard collision, there may be one or two or none (of flavor /) within the 
acceptance window. It is straightforward to obtain the first Ex-nioment for each flavor 
/ in the hard collision, 



0"hard(v / S,Po)(-Ex) = / dE T 



da f 
dE^ 

X 



/ dp\dy x dy 2 ^ Xifi/p&i, Q 2 ) x 2 f j/p (x 2 , Q 2 ) 1 
J 1 + Ok 

Pre(yi) (20) 



- + bki 

<*!> 



fk ~dl ^ ,U ' + fl ~dt ^ ' 
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i.e. basically as the integral of Eq. ([!]) over y in the accepted ^/-region, weighted by 
Pt- 

In an A A collision with an impact parameter b, the average transverse energy carried 
by semihard partons of flavor / to the acceptance region becomes then 

E^(b, y/s,p ) = T AA (b)a h; , ld (y/s,po)(EJ r ), (21) 

where T^(b)o"hard(\/SjPo) is the average number of semihard collisions (all y) and (E^) 
is the average transverse energy of the flavor / at \y\ < 0.5 in each semihard collision. 

3 The results 

The following results are based on a numerical evaluation of Eqs. (|]) and (pPf). For 
definiteness, we give the results for all y, in spite of the fact that very large \y\ phe- 
nomena cannot be experimentally observed in the collider mode and in spite of the 
fact that they are theoretically less well founded than the results near y — - they 
involve structure functions at extremely small x ~ p^/s, where they have not yet been 
measured, for nuclei, in particular. 

Firstly, to give an idea of the subprocesses contributing to g, q and q production, Fig. 
1 gives the y distributions from the various subprocesses at y/s = 5500 GeV. The scale 
at left refers to pp collisions, that on the right, obtained by using Eqs. ([12])- ([H]), gives 
the corresponding number of partons produced in an average central Pb+Pb collision. 
This scaling from pp to Pb+Pb is not valid for large \y\, where it may lead to a violation 
of energy or baryon number conservation. For kinematical reasons very small values 
of x are involved (see Fig. 4 below) and shadowing becomes important. In Fig. 2 
for \fs= 200 GeV one is far from violating energy or baryon number conservation in 
Pb+Pb for |y|<4 and one can thus apply the right hand scale at least up to this y; for 
pp there is no problem at any y. 

One observes that for gluons the process gg — > gg is dominant and increasingly so 
with increasing y/s. For quarks the dominant process is gq — > gq and the forward- 
backward peaking is due to valence quarks in this subprocess. For antiquarks the 
dominant process is gq — > gq, without any valence quark contribution. 

The y distributions are summarised in Fig. 3, which also shows the important net 
quark number distribution. Some numerical values are given in Tables 1 and 2. 
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Figure 1: The rapidity distributions of gluons, quarks and antiquarks (pp > Po — 2 
GeV) produced in pp collisions at ^/s = 5500 GeV together with the subprocess con- 
tributions. The axis to the right (to be applied only for |y|<4) gives the corresponding 
number of partons produced in an average central Pb+Pb collision. 

3.1 Initial state 

On the basis of the previous results one can immediately discuss the appearance of 
the system at the initial time 7* = 0.1 fm. For concreteness we give fixed numbers 
neglecting possible errors - and keep even too many decimals. Due to the expected 
dominance of semihard processes in the TeV energy range, let us mainly focus on the 
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Figure 2: As Fig. 1 but for yfe = 200 GeV 

results for LHC. In the figures and in the tables, we give the numbers for y/s = 200 
GeV for comparison, so that our discussion may be reproduced for RHIC as well. Thus, 
unless otherwise stated, the numbers apply to an average central Pb+Pb collision at 
y/s— 5500 GeV, p = 2 GeV and the rapidity range —0.5 < y < 0.5. The results in 
Fig. 1 and Tables 1 and 2 firstly imply that there are 

4350 gluons + 200 quarks + 190 antiquarks, (22) 

which carry a transverse energy of 

12960 GeV for gluons + 620 GeV for quarks + 590 GeV for antiquarks. (23) 
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Figure 3: Distributions of various quanta for y/s = 5500, 200 GeV. Note that the total 
hard and gluon cross sections have been scaled down by a factor 10 and 5 for yfs = 5500 
GeV and 200 GeV, respectively. 





range of y 


total 


9 


Q 


q 


5500 


\y\ < 0-5 


74.06 


67.95 


3.123 


2.988 




all y 


657.5 


555.8 


68.34 


33.06 


200 


\y\ < 0.5 


1.885 


1.556 


0.2025 


0.1266 




all y 


10.00 


7.077 


2.278 


0.6454 



Table 1: Values in mb of the perturbative parton cross sections Ohardlv^Po) an d 
^hardlv^'Pcb \y\ < 0.5) in pp collisions at y/s= 5500 GeV, p = 2 GeV. The number of 
partons of type / = g, q, q in an average Pb+Pb collision is Np hPh = 2crjJ ard Tp bPb (0) ~ 

64cT Ld- 



In a space-time picture these are formed at a time 1/po = 0.1 fm after the collision and 
thus, using eqs. QlBD , the corresponding number and energy densities are 

325 14.9 14.2 

%-^3> n<? -fm^' n? ~ta3' ( } 

GeV GeV GeV . s 

69 = 967 h^' e9 = 46 - 3 ftr^' e5 = 43 - 9 finT- (25) 
Note that here we consider an idealized (simplified) picture of a very high energy 
heavy ion collision, where the transit time of the colliding nuclei is reduced to zero. For 
the TeV energy range, this is a good approximation, since then 2i?p b /7 1/Po- How- 
ever, semihard and softer gluons extend longitudinally much further than the strongly 
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range of y 


total 


9 


Q 


Q 


5500 


\y\ < 0.5 


443 


405 


19.4 


18.4 




all y 


3830 


3227 


407 


196 


200 


\y\ < 0.5 


10.02 


8.211 


1.127 


0.6834 




all y 


51.99 


36.53 


12.09 


3.370 



Table 2: As Table 1 but for values of the first S T -moments cxhard 

{E{) in mbGeV of the 

parton distribution. The transverse energy of partons of type / in an average central 
Pb+Pb collision is = 0Wd(-Kr)^PbPb 

(0) «32(T hard (Eft. 



contracted valence components. For the semihard partons, the typical uncertainty in 
the longitudinal direction is Az ~ 2/ (xi( 2 )y/s) ~ 1/po- This width in the production 



mechanism |28| is also neglected in the idealized picture. Then, in the central rapidity 
unit, one may estimate (E) « 2(E^) sinh0.5 = 1.04(.&r) for the initial parton produc- 
tion, so that converting the transverse energy directly into energy density is a good 
approximation. 

Let us briefly recapitulate the discussion of the gluonic subsystem || , which has two 
notable properties. Firstly, the Et per gluon is 

^9 e pQCD 

^ = ^qcd =2.98GeV, (26) 
JV PbPb rig 

which naturally is somewhat larger than p . For an ideal gas of massless gluons in a 
complete thermal equilibrium, the temperature can be determined from 

£ ideal = 37l 2^ QQ . 16T< 4 ^ 

For an ideal gas with e 1 ^ 1 = e^ CD we find T eq = 1.10 GeV. It is interesting, and 
surprising, to see that 



^pQCD gideal 



pQCD ideal 
Tig Ug 



2.7T eq . (28) 



As far as energy per particle goes, the initial gluon subsystem thus is "thermalised" from 
the beginning. Note especially that we have not considered isotropisation of the system 
at all. In our idealized picture, (p 2 z ) (Pt)/2 initially in the central rapidity unit and 
collisions are needed to make the system isotropic. What our conclusion does indicate, 
however, is that since the appropriate energy/gluon is already there, no additional 
particle production is necessary, increasing thus possibilities for a rapid thermalisation. 
At this point it is perhaps worth mentioning that both the time of isotropisation and the 
conclusion of thermalisation depends on modelling of the initial collisions, especially on 
how the longitudinal width of initial parton production is accounted for. How rapidly 
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isotropisation then proceeds depends on the expansion of the system and on the rate 
of secondary collisions. Note that if there exists a finite width in the longitudinal 
direction, one may estimate the isotropisation time (and the initial conditions) also in 
a free streaming picture || 

Secondly, with the present p = 2 GeV the gluonic subsystem satisfies precisely the 
transverse saturation criterion (TB) for = 5500 GeV. In general, for pp collisions, 



Po is a phenomenological parameter separating the computable perturbative physics 
from the more model dependent soft physics. The value of Pq depends on what is 
assumed for the soft physics, and it is restricted by the measured total and inelastic 



cross sections [y| and transverse momentum spectra of charged particles in pp and 
pp collisions [|17]]. As we have shown above, sufficiently many semihard partons are 
produced in a heavy ion collision at very high energy, so that p gets a better defined 
dynamical significance, indicating that the soft component becomes less relevant. From 
Table 1 we can see that there is no transverse saturation at po = 2 GeV for a/s = 200 
GeV but the saturation occurs for smaller p , involving p^'s in the non-perturbative 
regime. In fact, saturation is reached at po = 1.6 GeV for y/s—1.8 TeV and at po = 1.0 
GeV for a/s=200 GeV. Therefore, at the RHIC energies, one should not neglect the soft 
component. A phenomenological way to obtain the saturation (dynamical screening) 
in a self-consistent manner in the initial parton production was suggested in Ref. |29j . 

Consider then the (anti)quarks, the density of which is below the thermal density by 
a factor of about 30. From Eq . (|24]) one sees that the initial net baryon number density 
at the time 1/po = 0.1 fm is 

1/ x 0.21 , , 

n B _s = -{n q - riq) = (29) 

Thus even at these ultrarelativistic energies the initial net baryon number density is 
more than the usual nuclear matter density, even if the initially baryon-rich fragmen- 
tation regions are already far apart in the forward and backward directions. Similarly, 
one may also compute the initial asymmetry: 

|— | = 18%.. ..29%, v^ = 200GeV, 
B + B 

= 4.7%. ..6.3%, = 1800 GeV, (30) 

= 2.1%.. .2.6%, = 5500 GeV, 

where the range of values shown corresponds to varying p from 1 to 3 GeV and 
where the result also is given for lower energies for comparison. Note that the initial 
asymmetry depends only mildly on p Q . The difference between yfs = 5500 GeV and 
200 GeV is clearly due to the increase in the quark and antiquark rates with increasing 
energy, the difference B — B varying less rapidly. Although the mechanism considered 
here is perturbative, the initial asymmetry for the LHC energies is of the same order 
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of magnitude as what is obtained in Ref. |23] in a purely non-perturbative gluon field 
model. 

Finally, since for a thermal boson gas s = 3.60n, the total initial entropy (using the 
number of effectively thermal gluons) is S = 15900 and, since B — B = 2.9 initially, 
the initial net baryon-to-entropy ratio at T, = 0.1 fm is 

^^1.8-10- 4 ~— . (31) 
S 5000 V ; 

It is also interesting to compare this number with the one for the early Universe, where 
the inverse of the specific entropy is ~ 10~ 9 . Thus we are still relatively far from the 
extreme conditions in the early Universe. 



3.2 Further evolution 

Further evolution of the system will go through possible thermalisation (quarks also 
may be thermalised) , plasma expansion, phase transition, hadron gas expansion and 
decoupling stages. It is clear that any conclusions will at present be model dependent. 

The simplest assumption is that of instant thermalisation and further essentially 
adiabatic expansion during all the stages, in analogy with the early universe. If the 
entropy (S ~ 15900), net baryon number (B — B m 3) and the net baryon number-to- 
entropy ratio ((B — B)/S = 1/5000) are all conserved, the initial values are also the 
final ones. For a hadron gas S"/iV>4 depends somewhat on the decoupling temperature. 
For a pion gas S ~ 4iV and a multiplicity of about 4000 is predicted. 

The assumption of instant thermalisation for quarks is, of course, unrealistic. In- 
stead, one may assume that they thermalise adiabatically at some later time r q . From 
the constancy of entropy it then follows that g g T 3 = (g g + g q )T 3 ew , i.e., 

T new = T(l + g q /g g y 1/3 ; (32) 

the increase g g — > g g + g q in the number of thermalised degrees of freedom is com- 
pensated for by a decrease in T. The total entropy and the multiplicity prediction is 
unchanged, but the energy density is decreased by the same factor as T. Note that 
the adiabatic scenario of quark thermalisation leads to a lower limit of the final pion 
multiplicity. Non-equilibrium models, where entropy, and therefore also the pion mul- 
tiplicity, increase due to the chemical equilibration of quarks have also been presented 

Assuming that each produced particle carries on the average a transverse energy of 
itlt ~ 0.5 GeV, the final transverse energy in the thermal scenario is about 2000 GeV. 
What happened to the large initial energy of 13000 GeV in Eq.([25|) in this scenario? 
It has been transmitted to larger rapidities by the work done against longitudinal 
expansion. In other words, the total entropy and number of particles in a unit of 
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rapidity remain constant, but since the energy/particle ~ T ~ 1/t 1 / 3 decreases, the 
total energy in a unit of rapidity decreases. The total reduction factor is the ratio 
between initial and final values of energy /particle = 2.7Ti/rriT ~ 6. 

We thus come to a series of scenarios starting from fully adiabatic expansion proceed- 
ing via increased amount of dissipation to a free-streaming expansion. The predictions 
of the final multiplicity and E-? starting from the same initial values in eqs. 
thus can vary from Na n = 4000, .Erfm = 2000 GeV to values larger than these by a 
factor of up to six. To narrow this range, reliable estimates for mean free paths would 
be needed. 

It may be of some interest to show how increasing dissipation can be described by 
the introduction of a shear viscosity |3(| f|. If one writes r\ = r] T 3 and uses e = 3aT 4 , 
the solution of the usual longitudinal hydrodynamical equations for T is 



T 



T, 



Ti 



1/3 + mJ_ 

t ) 3a 2rj 

1/3 



1/3 



Ti 
T 



6ar,- 



(33) 



Dissipation thus increases the temperature one reaches at a certain proper time and 
thus it also increases the lifetime and the multiplicity. All of this is parametrised by r/Q. 
By definition, the effect must be small, i] <C 6aTjTj. In the present context, r] could, 
in principle, be measured by observing deviations from the adiabatic multiplicity of 
about iV fin = 4000. In practice, however, this leading estimate is not reliable enough 
to permit measurements of corrections to it. 



4 Soft component? 

One of the motivations for the choice of taking po = 2 GeV was that at this value the 
produced gluons are so numerous that they are transversally saturated (Eq. flTBD ) so 
that []J this semihard component gives a quantitatively correct estimate of an average 
event at LHC energies. At lower energies the saturation is reached for still smaller 
values of p (at p = 1.6 GeV for ^=1.8 TeV, at p = 1.0 GeV for ^=200 GeV), 
for which next-to-leading order corrections are expected to be larger. At SPS energies, 
\fs= 20 GeV, the whole approach becomes useless. At these energies the bulk of the 
events is described by a nonperturbative soft component based on string or colour field 
formation and decay. 

A way to assess the importance of the corrections is to assume that the soft com- 



ponent, dominant at \fs= 20 GeV, is constant [16]. The SPS data then imply |T| 
that this soft component gives an average contribution of E^ ft (0) ~ 420 GeV in the 
central rapidity unit of central Pb-Pb collisions. The corresponding perturbative num- 
ber ^ ard (0) is 14 200 GeV at y/h= 5500 GeV (p = 2.0 GeV, Table 2). With these 
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numbers, (assuming that for the SPS the hydrodynamical stage of the system is not 
very long) a simple order-of-magnitude estimate for the final baryon-to-entropy ratio 
can be made. Let us decompose the ratio into soft and hard contributions: 

B-B (B-B\ (B-B\ 



S 



Phard 



Q—) + ftoft(— s~ , (34) 

O / hard \ D J soft 



where phard = S'hard/S' and p so ft = S so f t /S. An estimate for these is obtained from the 
transverse energies above, 



IP hard ipsoft 



Phard ~ p hard p so ft j Psoft ~ phard , psoft ' (^5) 



Using then the baryon-to-entropy ratio as in Eq. (|31| ) for the hard component, and 
1/50 from Ref. j[3] for the soft component, we obtain the final baryon-to-entropy ratio 

B — B a , 

— 8 • 10~ 4 (36) 

for nuclear collisions at LHC. 

For RHIC, the corresponding estimate is more unreliable due to the following reasons: 
Firstly, as discussed above, at \fs= 200 GeV, the saturation of minijet production 
occurs at px ~ 1-0 GeV but as this may already be in the region where the minijet 
cross sections become unreliable, let us consider a range from p = 1.5 GeV to 2.0 
GeV. This gives E^ ld (0) = 708 GeV to 326 GeV. Secondly, although the gluons do 
dominate the early parton system also at RHIC, the energy /gluon is not thermal, so 
the assumption of instant thermalisation of glue at r ~ 0.1 fm is probably not as good 
as for the LHC. Following, nevertheless, the same line of reasoning as for the LHC 
gives (B - B)/S ~ 2 • 10~ 3 initially at r ~ 0.1 fm, and (B - B)/S ~ 9 • 10~ 3 for the 
final baryon-to-entropy ratio for nuclear collisions at RHIC. 

We see that the nonperturbative sources may increase the baryon-to-entropy ratio 
considerably, by a factor of about 4, from the initial values both at LHC and RHIC. 
At the LHC, the increase is dominantly due to the increase in the net baryon number 
and not in the total entropy, while at the RHIC both the net baryon number and the 
total entropy are increased by the soft source contributions. 



5 Very forward and backward partons 

The previous discussion has been confined to the y = region, although the curves 
in Figs. 1-3 are computed over the entire rapidity range. They are reliable on the 
pp level, but the simple transformation in Eqs. ( |i~2| , |i"4"D used to transform them to 
the A+A level assumes independent subcollisions and thus does not take into account 
energy and baryon number conservation. It thus fails at large y, where a few gluonic 
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sub collisions can consume all the available energy or valence quark+gluon sub collisions 
can consume all the baryon number. For the approach to be valid one has to be far 
from these limits; there is no problem at RHIC, but at LHC one should not go beyond 

\y\ ~ 4. 

Another problem is that kinematically a very forward parton demands that the 
other one also be at large y. This, on the other hand, implies that the x value of one 
of the partons is extremely small and in a region in which the structure functions are 
unmeasured. This effect is quantitatively discussed in Fig. 4. 

6 Conclusions 

We have in this paper computed the initial number of gluons, quarks and antiquarks 
with p T > p = 2 GeV (value determined by transverse saturation) produced in Pb+Pb 
collisions at y/s= 5500 GeV near y = and studied the further evolution of this system 
especially in view of the behavior of the net baryon number. The main conclusions 
are that the produced gluon system is transversally saturated and the parameter p 
gets a dynamical significance in LHC nuclear collisions. The perturbatively produced 
gluon system has also energy/gluon as in a thermal ideal gas, so only isotropisation 
is needed to make the system thermalised. This clearly enhances possibilities of a 
very rapid thermalisation. The initial net baryon number density at Tj = 0.1 fm may 
be of the order of nuclear matter density, or even larger, but this is rapidly diluted. 
Initially at r { = 0.1 fm the baryon-to-entropy ratio is (B — B)/S ~ 1/5000, and the 
final baryon-to-entropy ratio may be larger by up to a factor 4. We conclude that 
the baryon-to-entropy ratio is small, though still orders of magnitude larger than the 
corresponding quantity in cosmology. 

The main sources of uncertainly in the computation are next-to-leading order cor- 
rections, nuclear shadowing of the initial parton distributions and dissipative effects 
and the non-perturbative source terms during expansion and hadronisation. Next- 
to-leading order corrections with the acceptance conditions applicable here have not 
been computed but should be. Nuclear gluon shadowing can only be modelled and 
ultimately has to be measured. This also holds for dissipative effects. Although the 
detailed numbers are uncertain, we believe that this general framework will be used to 
analyse the data when it finally comes. 
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Figure 4: The forward-backward peaking of the quark distribution in Fig. 1 
arises from valence-gluon subcollisions. This figure shows the parton luminosities 
Xig(xi, p )x2V(x2, Po) and xiV(xi, Po)x2g{x 2 ,Po) as a function of the rapidity y 2 of 
the gluon, which is integrated over in Eq. ([Tp. The X\ and x 2 values of the initial va- 
lence quark or gluon are also shown as a function of the rapidity y 2 ■ For very forward 
quark rapidities (y = 6) the dominant configurations are such in which also the final 
gluon is at y 2 ~ 6, but then also the x 2 is extremely small (~ Po/s) and in a region 
where the structure functions are unknown. 
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